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I .  IMTRODUCTION 


This  report  is  Intended  to  serve  three  purposes.  The  first  Is  to  provide 
soa«  discussion  of  the  structure  and  functioning  of  the  niolecular'*%lectromag'* 
netlc  correlation  foraallsa  (MECF)  which  Is  applicable  to  the  study  of  gases 
Interacting  with  the  electromagnetic  field.  The  second  Is  to  present  an  anal** 
ysls  performed,  by  the  author,  during  77  1986  In  the  application  of  this  for- 
mallsm  to  the  study  of  pressure  broadening  In  gases.  The  third  is  to  briefly 
summarize  past  accomplishments  and  some  possible  future  applications  of  the 
approach.  It  is  hoped  that  the  reader  will  gain  some  understanding  of  the 
nature  and  scope  of  this  rigorous  and  comprehensive  approach  to  the  descrip*- 
tlon  of  gas/electromagnetlc-fleld  Interactions^. 

Because  application  of  the  general  theory  to  pressure  broadening  in  gases 
Is  only  In  Its  Initial  development  state,  the  analyses  presented  in  Section 
III  Is  only  one  milestone  in  the  study  of  the  Inclusion  of  Intermolecular 
interactions  due  to  Coulomb  forces  in  the  physical  description.  The  analysis 
presented  is  qualitative  in  nature,  and  develops  the  techniques  and  data 
needed  to  evaluate  the  perturbations  Induced  by  increasing  the  density  and/or 
pressure  of  the  gas  in  the  sample  volume  considered.  Because  of  its  simpllc~ 
Ity,  the  hydrogen  atom  has  been  chosen  as  the  "molecular”  species  used  in  this 
initial  Investigation.  It  Is  hoped  that  many  of  the  qualitative  and  semlquan* 
tltatlve  features  of  the  phenomena  of  Interest  can  be  determined  by  the  use  of 
this  "test  molecule." 

Finally,  as  a  brief  Indication  of  the  comprehensiveness  of  the  MECF  for* 
mallsm.  Section  IV  outlines  several  past  applications  of  the  method  and  points 
out  some  potential  areas  for  future  applications;  these  latter  areas  range 
from  a  study  of  the  dynamics  of  the  free^electron  laser  to  consideration  of 
adverse*eondltlon  Imaging  in  fog  or  mist.  It  Is  hoped  that  this  discussion 
will  stimulate  Interest  In  the  capebllitles  of  the  MECF  approach. 


^C.  A.  Coulter  and  D.  W.  Howgate,  Physical  Review  A  31,  843  (1985). 


II.  BACKGROUND 


Th«  MECF  formalism  is  a  comprehensive  and  rigorous  method  for  describing 
the  interaction  of  a  molecular  gas  with  the  electromagnetic  field.  It  con¬ 
tains  only  the  very  mild  physical  approximations  (see  paragraph  A  below).  Up 
to  this  point  the  MECF  description  has  been  used  in  a  form  that  applies  to  a 
sunosMlecular  gas,  but  extension  to  a  mixed  gas  is  possible.  The  derivation 
of  the  MECF  description  requires  the  use  of  quantum  field  theory.  A  key  step 
in  the  development  of  the  theory  is  the  introduction  of  an  operator 
that  can  be  thought  of  as  the  annihilation  operator  for  a  molecule  with  inter¬ 
nal  state  k  and  center  of  isass  x_.  However,  the  MECF  formalism  itself  is 
expressed  solely  in  terms  of  non-operator  expressions.  It  contains  the  energy 
eigenvalues  for  the  molecular  species  of  interest  and  also  certain  functions 
that  are  calculated  from  the  wave  functions  for  the  molecule.  The  most  impor¬ 
tant  of  the  functions  determined  by  the  wave  functions  are  those  shown  in 
paragraph  B  below.  The  most  important  of  the  dependent  variables  are 
discussed  in  paragraph  C. 

A.  Physical  Approximations 

Physical  approximations  contained  in  the  MECF  description  of  the 
electromagnetic  interactions  of  a  molecular  gas  are: 

•  Nuclei  are  treated  as  elementary  particles. 

•  Ihiclei  and  electrons  are  treated  nonrelativistically. 

•  The  effects  of  molecular  collisions  with  the  walls  of  the  container 
on  the  internal  anlecular  states  are  neglected. 

B.  Molecular  Wave  Functions 

The  MSt  important  functions,  appearing  in  the  MECF  formalism,  are 
determined  from  the  laolecular  wave  functions.  Definitions  of  these  functions 
are  given  in  the  following  expressions.  In  these  expressions  is  a 

function  that  specifies  the  most  important  part  of  the  coupling  of  the  mole¬ 
cules  to  the  electromagnetic  field,  •{(•(jc)  is  a  function  that  describes 
the  intermolecular  interaction,  z  is  used  to  represent  the  complete  set  of 
internal  coordinates  for  the  molecule,  and  Uf^lz)  is  the  wave  function  for  the 
internal  molecular  state  with  index  k. 

Sjti(£)  -  ^  {ftajiqi/2i 

ij 


and 

Vkll'k’^i)  <li<li' 

ii' 


•i)  J  5(£-;rj)(u*(z)7jUj(z) 

-  (7ju*(z))u  (z)lTTd3zj 

VC  rv 

J 

J"  (u*(z)u*(z’)Uj,(r')U|^,(z)/|  x^-(2^-^')  j  ) 

J  y 


(1) 


(2) 
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Here  x  la  the  oolecular  center  of  ouiss,  displacement  of  the  1-th 

particle  from  the  center  of  mass,  a^^j  Is  the  matrix  that  transforms  laboratory 
coordinates  to  Internal  molecular  coordinates,  and  and  mj^  are  the  mass  and 
charge,  respectively,  of  the  1-th  particle  In  the  molecule. 

C.  Dependent  Variables 

The  dependent  variables  In  the  MECF  formalism  are  a  set  of  correlation 
functions  for  the  molecules  and  for  the  electromagnetic  field.  The  most 
Important  of  these  are  the  low-order  MECFs  expressed  as  expectation  values  of 
the  electromagnetic  vector  potential  A(x,t)  and  the  molecular  annihilation 
operator  adjoint.  V  Is  the  volume  of  the  gas,  and  C}}  is  a 

constant  dependent  on  the  number  N  of  molecules  in  the  volume  and  on  the 
number  of  particles  of  each  species  In  a  molecule.  An  expression  for  the 
value  of  Is  given  In  the  paper  referenced  in  Section  I. 


A(x,t)  -  <A(x,t)>,  P^j^CXj^.x^.t)  -  VCjj<P+(x^,t)Pj^(x2,t)>, 


F«b(x^,X2,t)  -  <Aa(x^,t)Ab(x2,t)>,  H«b(Xj^,X2,t)  -  c'2<  j^A^Cx^ ,t)  3^Al>(x2,t)>, 


®***^il»i2**^^  *  (2c)‘*kA«(Xj^,t)3^A^(x2,t)  +  3^A^(x2,t)A*(Xj^,t)>, 


t*3»54»  t)  -  VCjj<Aa(x^,t)Ab(x2,t)P+(x3,t)Pj(x^,t)>.  (3) 


Calculation  of  the  functions  and  shown  In  paragraph  B  above 

Is  a  sizeable  task  even  for  the  case  of  the  hydrogen  atom,  where  the 
"molecular"  wave  functions  are  available  In  analytic  form.  Determination  of 
the  functions  in  other  cases  will  be  even  laore  demanding,  but  Is  not  out  of 
the  question  (at  least  in  approximate  form)  for  some  of  the  other  simplest 
atoms  and  molecules.  Figure  1  gives  the  status  of  the  "data  bank"  of  these 
functions  as  of  October  1985.  The  functions  T  and  U  referred  to  In  this 
figure  are  required  for  the  description  of  much  weaker  Interactions  of  the 
molecules  with  the  electromagnetic  field  (Involving  polarization  and 
effects),  and  have  not  yet  received  much  study.  ~ 


MECF  DATA  SANK  (E,  S.  T.  U.  V)  ~  ^ 
FY  85  -  87 

Targeted  Species 

e  -  electron 
H  -  hydrogen  atom 
H2  -  hydrogen  molecule 
He  -  helium 


SPECIES/PARAMETER  MATRIX 


Definition  of  Symbols  Utilized  (not  unique  to  MECF) 

0  -  no  requirement 

F  -  requirement  trivial. or  available  from  standard  sources 
I(N)  -  requirement,  but  no  input  as  of  end  of  fiscal  year  198N 
P(N)  -  requirement,  with  partial  input  as  of  end  of  fiscal  year  198N 
F(N)  -  requirement,  with  full  input  as  of  end  of  fiscal  year  198N 
E  -  energy  of  targeted  species 


Figure  1 .  Sui—ary  of  the  MECF  data  bank  as  of  10  October  1985. 
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MECF  Parameter  Matrix  Elements 

S  -  principal  field/matter  parameter 
V  -  matter/matter  parameter 
T,LI  -  apply  only  at  very  short  wave 
lengths  and  very  high  fields 


MAHER /MATTER  MATRIX  (V) 


e  H  H2  He 


The  MECF  formalism  reduces  to  the  standard  semlclassical  approximation,  if 
one  assumes  photon  wavelengths  are  significantly  greater  than  molecular  dimen¬ 
sions  and  low-molecular  number  densities  so  that  the  intermolecular  interac¬ 
tion  terms  can  be  neglected.  In  addition,  it  is  necessary  to  assume  that  the 
quantum  correlations  of  the  electromagnetic  field  with  the  molecules  can  be 
neglected,  so  that  approximately 


4i<i’  Xj,£2*t)  *  A(x.»t)  ^kl 


(4) 


With  these  approximations,  one  is  led  to  a  closed  set  of  equations  involving 
only  and  the  functions  and  internal  molecular  energies 

At  higher  gas  pressures,  the  effects  of  the  intermolecular  interaction 
term  involving  the  functions  V|(2nin  must  be  included.  The  corresponding  set  of 
equations  also  contains  the  new  MECF  Pklmn*  and  to  avoid  having  to  extend  the 
set  of  MECF  equations  significantly,  it  is  necessary  to  approximate  this  four- 
point  correlation  function  in  terms  of  the  two-point  functions  Any  such 

approximation  must  take  into  account  the  properties 

S/Pkk(  £,3t,t)d3x  ■  V, 
and 


where  the  or  -  sign  is  chosen  in  the  last  equation  for  an  even  or  odd 
number  of  fermions  per  molecule,  respectively.  An  approximation  for  Pkimn 
that  nearly  satisfies  the  first  and  third  of  these  conditions  (with  an  error 
in  satisfying  the  first  condition  of  relative  size  equal  to  the  ratio  of  a 
molecular  volume  to  the  total  volume  of  the  gas)  is 


(7) 
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This  approximation  is  employed  in  the  following  analysis  of  intermolecular 
interaction  effects  to  avoid  the  introduction  of  higher  order  MECFs  at  this 
point  in  the  development  of  the  implications  of  the  theory. 


Note  that  the  quantities  necessary  to  describe  the  physical  properties  of 
the  system  can  be  expressed  in  terms  of  the  MECFs.  For  example,  the  expecta¬ 
tion  values  of  the  molecular  number  density,  molecular  momentum  density, 
Poynting  vector,  and  electric  and  magnetic  energy  densities  can  be  written  as 


n(x,t)  -  (N/V)XPkk(x,x,t), 


£(x,t)  -  (N/V)2^<‘/2i)((7x-7xi)Pkit(x*,x,t)I  (  x'-x’ 


<S*(£,t)>  -  cJ3(7^G**>(x,£',t)  -  7«G*»*’(3t,£',t)]  (  x'«x’  (8) 


<l/2(E(x,t)]2>  -  l/2l3H**(x,x,t), 
—  —  a  —  “ 


<l/2[B(x,t)]2>  .  1/2  Z  Fa*(x,x’,t>  -  7^7*  Fab(x,x ' , t) ]  |  x'-x. 


In  deriving  and  illustrating  the  properties  of  the  MECF  formalism,  it  is 
desirable  to  use  a  simple  "molecular"  system  for  which  the  calculation  of  the 
necessary  energy  values  and  wavefunction  dependent  functions  ^i  is  as  simple 
as  possible.  The  hydrogen  atom  is,  of  course,  ideally  suited  to  this  task. 
Therefore,  it  is  used  as  the  illustrative  system  in  the  succeeding  analysis  of 
intermolecular  interaction  effects. 


III.  ANALYSIS:  PRESSURE  BROADENING 


The  information  and  physical  properties  sought  in  the  implementation  of 
this  discussion  can  be  directly  correlated  to  an  analysis  of  the  current  vec¬ 
tor  J(xyt)  that  occurs  on  the  right  side  of  the  wave  equation 


(v2  -  c-232)  A(x,t)  -  -(kc)‘kN/V) 


(9) 


where  N/V  is  the  matter  density  parameter  and  the  initial  applied  electromag¬ 
netic  field  is  assumed  to  have  a  plane  wave  representation  of  single  wave 
vector  k  >  kic.  Since  utilization  of  the  hydrogen  atom  was  chosen  as  the 
"molecule,"  the  three  transitions  can  now  be  chosen  that  will  produce  the 
Lyman-a  line  as  the  input.  In  the  MECF  formalism,  these  three  transitions 
produce  field/matter  coefficient  vectors  (x)  that  lie  along  the  z,  x,  and  y 
axes,  respectively.  However,  if  the  direction  of  propagation  of  the  initial 
field  k  is  chosen  to  be  the  z-direction  and  its  vector  amplitude  A  to  lie 
along  the  x-direction,  then  J  will  also  lie  along  this  direction  and  now  the 

only  concern  is  with  that  upper  level  (labeled  2)  which  supports  an  sj  along 

“Im 

the  x-direction.  In  this  reduced  space  one  can  now  express  J(x,t)  as 


J(x,t)  -  2k3c  Re(X(t)A(x,t)I,  X(t)  •  (fiA)'kkc)'2  s(jc)2  (0,t),  (10) 

where 

S(k)  ■  i'l^d^xexpC  i_k*x)  “  -i(hAc)”l 

S(k)  A(R,t)  P®2(£,t).  (11) 


In  arriving  at  this  representation,  expression  of  the  principal  molecular 
correlation  function  as  P;(]^(^,r,t)  was  chosen  in  terms  of  the  coordinates 
^(2Cj+X2^^2»  equation  of  motion  for  Pj2^— been  effec¬ 

tively  linearized,  by  assuming  that  the  diagonal  elements  PrmnC  t)  which 
appear  on  the  right  hand  side  can  be  approximated  by  time  independent  repre¬ 
sentations  gn|P°(£)  which  are  also  independent  of  the  global  variable  R.  The 
parameters  a  and  A  remain  to  be  chosen,  whereas  the  function  X(t)  is  assumed 
to  exhibit  the  form  of  a  volume  susceptibility* 


Operating  from  the  left,  on  the  off-diagonal  element 

0'(R,r,t)  -  i3t  -  «i/M)VRi7r  -  (E'2-E' (12) 
and  only  the  principal  field/matter  interaction  terms 

i(ch)-ls(k)  [gjA(R-r/2,t)  -  g2A(R+r/2,t) ]P«(£) , 

A(3c,t)  “A  exp( i(^.x^-kct)  1  (13) 

with  the  Coulomb  matter/matter  terms 

(^)(N/V)f"lJ*d3x  ^2122^—^  (82^22^—”—^^’—^—*^^ 

■  82^22^— 

that  appear  in  the  MECF  formalism  will  be  evaluated  on  the  right-hand  side  of 
the  ^22^— equation  of  motion.  The  energies  are  represented  as  in 

order  to  imply  that  they  may  contain  additional  elements  brought  over  from  the 
right-hand  side  of  the  equation.  The  matter/matter  interaction  elements  have 
the  following  syimetry  properties:  ''^imni^i^  *  '^jmni^i^  “  '^inmj^i^  *  Ymijn^"*.^* 
The  intent  here  is  to  emphasize  the  qualitative  nature  of  the  analysis  and  to 
minimize  detail  whenever  possible.  With  this  in  mind,  the  consideration  of 
terms  on  the  right-hand  side  containing  only  diagonal  elements,  which  are 
independent  of  the  existence  of  an  applied  field  and  introduce  corrections  to 
Pniin(  R,£,  t)  of  order  (N/V)^,  is  left  to  the  Appendix.  The  above  treatment  is 
consistent  with  general  intent  to  depict  the  impact  of  the  Coulomb  interaction 
as  an  expansion  in  powers  of  the  molecular  density  parameter  (N/V)  or  its 
equivalent. 

Now  the  parameters  a  and  A  can  be  developed  according  to  the 
specifications 

P«(£)  -  exp [-(or) 2/21,  a  -  (MkBT)l/2/^^ 

Vii22(i)  -<iVA'/*d3k  V2222<ii)»  ^  -  (2ir2fiV)”lN(q  a  )2,  a'  -  A/N,  (15) 

where  T  is  the  temperature,  M  is  the  hydrogen  atom  mass,  and  q^  and  a^  are  the 
electron  charge  and  first  Bohr  radius,  respectively.  In  so  doing,  72222^*^) 
becomes  a  dimensionless  quantity  and  the  parameters  a  and  A  have  the  dimen¬ 
sions  of  inverse  length  and  inverse  time,  respectively.  The  choice  of  a 
simple  Gaussian  for  P®(£)  is  a  natural  consequence  of  the  requirements 
for  the  noninteracting  gas  to  be  initially  (t-0)  in  thermal  equilibrium.  It 
further  satisfies  the  additional  conditions  relevant  to  including  a  Coulomb 
interaction  under  a  perturbation  expansion  in  powers  of  A  or  N  (see  Appendix). 
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Nov,  for  convenience,  formally  remove  a  by  transforming  to  a  reduced  dimen¬ 
sionless  representation  by  way  of  the  relations 


♦  r ,  ff(MkBT) "1/2^  ■►  £,  (kB)T/M)  +  t, 


(kBT/M)“l/2k-l6  >  6, 

(kBT/M)“l/2ii-lA  ♦  Hr),  Pl2^1*t>*  1*  tk,  _6-6k,  (16) 

where  8_  is  an  arbitrary  wave  vector  and  one  is  now  normalized  relative  to  the 
influence  of  thermal  (kBT)  and  impulse  (hk^/M)  energy.  It  is  this  represen¬ 
tation,  by  way  of  Pi2(^»t),  that  will  provide  detailed  information  of  such 
physical  properties  as  phase  and  line  shape. 

Looking  at  the  equation  of  motion  for  o*'®  sees  that 


0(£,t)  Pi2^£,»t)  ■  Ri2^£.»t)»  0(£,t)  ■  6/A  +  iA"k3t  “  ^’^r)  •  (17) 

I  where  R22(£»t)  is  what  remains  on  the  right-hand  side  of  the  equation  after 

j  the  dependence  has  been  factored  out,  and  6  is  just  the  initial  detuning 

parameter  that  arises  from  expression  (12)  when  taking  the  derivative  of 
k(R,t)  with  respect  to  time  and  subtracting  the  energy  level  difference.  The 
inverse  of  0(^, t)  can  be  performed  through  the  use  of  the  appropriate  Green's 
I  function  relations 


Pl2(£,t) 


D(£, t-t')  Rl2(£.»t'), 


D(x,t) 


®i2^iL"£'  ;t), 


(18) 


a  representation  in  which  D(x^,t)  is  introduced  first,  leaving  the  final 
integration(s)  until  later.  To  properly  exploit  D(3C, t),  it  is  convenient  at 
this  time  to  introduce  the  function 


Q(£;^)  ■  8xpl-b2/2  ♦  U.Jb],  Q(£+i3c;£)  Q(^;x)  -  Q(£;0  Q(^+i£;£),  (19) 

where  £  and  £  are  two  arbitrary  vectors,  along  with  its  appropriate  integral 
representation  (transform) 


U(£,£)  ■  Tt"2/2y-l^j3,jl  Q(£;^-X»). 


(20) 
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The  following  loportant  relations  governing  the  application  of  both  D(x,t)  and 
U(8^,x)  can  now  be  developed:  ~ 

D(x,  t )  Q(  pa ;  u“lx)  -  Q(  u( 8+« ;  p-^t )  Q(  ps+i  p“lt ;  p'^x) 

-  exp(-l^*x^)  Q(p(8+^);p“l(t4»)) . 

UCpa.p-lx)  D(x,t)  Q(pa';p**lx)  •  V(p)“l  exp[-l^*x-l/2( y '8")2j  (21) 

•Q(p’(8++_6);(p')“l(t+x)), 

D(x,t)  UCps^jP^lx)  D(x,t')  Q(p£';p‘*^x)  -  V(p)“l  expt-l^*x;-l/2(  p’s")^) 

•Q( p • ( 8++^) ; ( p ’ ) “1 ( t+t ' +x ) ) , 

where  si*  (a,+  £')/2,  and  p'  •  2^/2jj  g^  arbitrary  scaling  parameter. 

After  inspecting  equations  (18),  (19),  and  (21),  it  was  concluded  that  as 
a  correlation  function  Pi2(£,»t)  is  important,  primarily,  within  a  radius  r<l. 
Also,  that  Pi2(£»t)  and  Ri2(£,. t)  change  rapidly  during  a  brief  initial  tran¬ 
sient  period  t<l,  and  then  approach  values  Pi2(£)  and  Rl2(£)  that  are  basi¬ 
cally  independent  of  time.  Thus,  for  times  t»r,  which  are  of  primary 
concern,  one  can  trrlte 


Pl2(l)  -  -lA  fdt  D(r,t)  Ri2(£).  (22) 

Jo 

where  Ri2(£)  nan  be  expressed  simply  in  dimensionless  representation  as 

M2(£)  -  -[giQ(-k/2;r)  -  g2Q(y2;rl 

(+)lx3/2vjd3k’  [gjU(k'+k/2,r)  -  g2U(k'-k/2,r) J 

•A  J  dt  D(£,t)  Ri2(£).  (23) 

with  the  aid  of  the  expressions  in  (21).  Using  relations  In  (21),  (22),  and 

(23),  it  is  now  easy  to  write  Pi2(£)  as 

Pl2(£)  ”  lAexp(-l^*r)J^dt  [g^Q(-y2+^;t+r)  -  g2Q(k/2+^;  rfr)  ] 

(+)x3/2^2g^(»j^.pj J*d3k'  ^^^22^^’^^^**^ 

•[g2exp(-l/A(k-tk’)2)Q(2l/2(^./2+^);2“l/2(t-H.+r))  (24) 

-  gjg2exp(-l/4k'^)Q(2^/2((^f4^)/2+^);2‘’l/2(^^^f+^))  J  ^ 

where  two  additional  terms  Involving  g.g2  and  g2  have  been  omitted  to  minimize 
the  amount  of  detail  displayed.  2 
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It  is  proper  here  to  think  of  Q(a^;M  as  one  of  the  primary  building  blocks 
developed  by  this  representation  of  the  general  formalism.  As  a  product  of  an 
oscillatory  element  exp(ia.b)  and  a  Gaussian  damping  term  exp(-b^/2),  its 
integration  over  time  allows  the  effect  of  damping  to  smooth  out  the  general 
6*character  produced  by  the  direct  time  averaging  of  oscillatory  motion. 
Furthermore,  as  far  as  the  important  Coulomb  interaction  element  Vj^22 
series  expansion  terms  of  arbitrary  order  are  concerned,  Q  (apart  from  a 

simple  multiplier  exp(-i_5.r))  uniquely  displays  the  functional  variation  of 
both  r^  and  t.  For  examp leT  for  terma  of  order  in  the  Vjj22  expansion, 

where  the  operator  U  has  been  applied  n-l  times,  one  finds  Q  displayed  as 
follows : 

Yn(_Q,jr)  exp( •f dtj.. 

o 

where  £^*£2 ■*■•••■*■£„ «  ii  arbitrary  vector  involving  and 

elements  of  both  Jc  and  the  Coulomb  interaction  wave  vector  1^'.  This  display 
is  a  transformation  from  the  basic  Q  representation  into  a  format  %rhere  a 
single  Gaussian  is  weighted  by  the  defined  function  Y,^(jl,0 .  This  function 
can,  in  turn,  be  displayed  by  using  the  complex  plane  to  unwind  the  multiple 
integration  over  time,  as 

n-l 

YnCn.O)  -  2  Yn-m(O.O) 

m“0 

^ds^  . . . 

where  Rn(^)  >-•  the  remainder  function  for  the  truncated  series  expansion  and 
introduced  in  expressions  (25),  (26),  and  (27),  is  just  the  projection  -'f 
the  vector  ^  along  the  direction  of  propagation  k.  The  elements  Yn(0,0)  form 
a  series  of  converging  constant  terms  Yj*(w/2)^/^,  Y2*l,  etc.  As  correlation 
is  of  secondary  interest  only  in  this  discussion,  r  has  been  set  to  equal  0  in 
developing  expressions  (26)  and  (27),  for  the  sake  of  simplicity  of  represen¬ 
tation.  For  Q^<1,  R(n)()can  be  expressed  in  a  Maclaurin  series  expansion  as 


W 


♦  (i"/n!)(u2n^)"Rn(tJ^), 
-n-l 

y*  ds"  exp((  un*”)^/2) , 


m 

fdtn  Q(urxn;(un>"^(£^-M>). 


Rn(£lk)  ■  1  >5^  (uj)"(2(m-l)*l)...(2*0^1)(n!/(n*2m)!)(n|^)2"',  (28) 

m«l 


providing  a  series  representation  of  infinite  extent  for  Y„(Q,0).  One  can 
also  define  the  generalised  moments 


M„P  -y*dQ^(Q^)P„p(.y^2n^2/2)  Yn(n,0)  -  ( -i )’’ J*dt  1 . . . 


/*dtndP(t*)exp(-(t*/iin)^/2), 
n  n 


where  6**(s)«(d^/ds**)  6(  z) ,  to  provide  an  additional  tool  for  analysis. 
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Th«  final  puah  comma  whan  ona  aceaapta  to  Incorporata  axpraaalona  such  as 
(25)  (or  ehoaa  of  ordar  n)  Into  thalr  total  snvlronaant,  which  nay.  In  turn, 
Inwolwa  additional  products  of  Gaussian  alaaants  containing  k,  and 
k**  ^  ganaratad  by  axprassions  such  as  thoaa  shown  liT  (Tl),  and  tha 

Couloab  iMighting  function 


(30) 


which  now  introducaa  tha  lapact  of  tha  particular  aolaeula  Involvad.  An 
attaapt  was  aada  to  structura  this  analysis  so  that  this  task  was  put  off  as 
long  as  possibla  >  thua,  optialslng  tha  ganarallty  of  tha  approach.  Although 
asprasslon  (24)  dapicts  Q  and  Its  ganaral  anvlronaant  at  first  and  sacond 
ordar  in  a  only,  it  is  basically  illustratlva  of  tha  principal  action  and 
procasaaa  that  ara  involvad.  Satting  r_ ■  0  and  insarting  ralatlon  (25), 
axprassion  (24)  can  ba  rastatad  slaply  as 


Pl2(0)  -  lAlgiTi(fr*/2,0)a*p(-(6-k/2)2/2)  -  g2Yi(^/2,0)axp(-(  5^/2)^/2J 
(>)w3/242y’d3iit  Vii22<Jl’)  •«l>I-lMk^^l/2(kJ-*-«Ht/2)2] 

•  t  g2T2  (i^'  /  2 ,0)axp(  -(  «-k/2  )  2  ) 


(31) 


-  *«i/2)V2)) , 


an  axpraasion  whara  tha  (4-)  tarns  display  tha  Gauaslan  product  undar  tha 
influanca  of  tha  Coulonb  Tntaraction  walghtad  by  tha  product  Y]^|^22^2 
avaragad  ovar  Ic' .  Tha  wava  vactor  It,  whata  It  appaars  as  an  lapluse  corrac- 
tion  dua  to  photon  anlsslon  or  absorption  (as  in  tha  first  tarn  of  (31)), 
should  hava  llttla  affact  on  tha  iapact  of  prassura  broadaning  and  rslated 
phanonana  at  wavalangths  far  axcaading  tha  Bohr  radius  of  tha  hydrogan  atoa 
and  at  finita  tanparaturas.  It  could,  in  fact,  ba  onlttad  fron  expression 
(31)  as  a  further  siapllf ication.  Tha  variables  k'  and  k'  that  appear  in 
expression  (31)  refer  to  tha  conponants  of  tha  vactor  It'  along  the  k  (or  z) 
direction  and  on  tha  x>y  plana  that  is  oomal  to  it,  respect ivaly. 


This  discussion  ends  by  ignoring  tha  Initial  population  in  the  upper 
level  (set  g2  *  0)  and  any  sinpla  iapulsa  corrections  (sat  Ic  >  0,  except 
where  it  occurs  inside  of  6).  Now,  tha^ratio  1(6)  can  ba  fomad,  using 


axpraasion  (31),  of  tha  tarn  of  ordar  6^  to  tha  tarn  of  order  6,  such  that 


1(6)  -  (?)1(7ii/9)a(H-1(2/i.)1-'2  sRi(  6))”^  /  d3k'exp[-l/4(k’2+k'2)-i/2{  5+k')2] 

J  *  y  r 


.2w.-2ji+(2/3)2(«*k')2j-«(l+l(2»)l/2( 5+k'/2) 

s  z 


•k’^k 

X 


-2(6+kV2)2R2(6-HcV2) 


(32) 
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wh«r«  k^,  ky,  and  k'  arc  tha  coaponanta  of  k ' ,  along  tha  x,  y,  and  z  axes, 
raapacclvaly,  and  a^  la  cha  cransforaad  (dltaanslonlass)  value  of  the  hydrogen 
Bohr  radlua  a^.  Aa  a  Uniting  approximation,  tha  renalnder  function 
will  satisfy  tha  relation 

Rti(z)  <  n2‘"/2zl/2/^b8(z).  (33) 


Using  0.121  nlcronatars  as  tha  wavelength  of  tha  Lynan^a  Una,  one  Is  now  In 
tha  position  to  evaluate  tha  aaqiUtuda  scaling  factor  &  as  a  function  of 
pressure  and  tanparatura.  With  tha  aid  of  expressions  (IS)  and  (16),  ^  may  be 
written  In  tha  fom 


A  -  3,80x10*  P/t3/2,  (34) 


where  pressure  P  la  In  atnoepherea  and  temperature  T  is  in  degrees  Kelvin.  At 
standard  condition  (P*l,  T*300),  then  A  “  0.732  a  value  slightly  below  one. 
Thus,  A  appears  at  first  glance  favorable  as  an  Initial  choice  of  a  scaling 
parameter  and  to  support  the  expansion  of  the  Coulomb  Interaction  as  a  pertur¬ 
bation  aeries  Involving  this  parameter.  This  completes  the  Initial  analysis 
where  certain  of  the  (t)  terms  arising  out  of  the  Coulomb  Interaction  have 
been  developed  within  The  confines  of  the  general  formalism.  This  analysis 
Indicates  that  proceduras  can  be  developed  for  measuring  the  effectiveness  of 
the  general  formaUsm  In  maetlng  moat  requlreamnta  In  physical  applications 
such  as  pressure  broadening,  etc.  The  further  development  of  such  procedures 
will  be  a  basic  part  of  the  next  phase  of  this  study. 


IV.  OTREI  APPLICATIONS  AND  ACCOHPLISHHENTS 


Kt  «n  integral  part  of  prior  davalopaant  of  the  MECF  approach,  the 
following  applicationa  have  bean  briefly  explored: 

1.  At  the  lowest  (aaaiclaaaical/low  density)  level  of  approximation,  the 
impact  of  gas  dynamics  on  the  stimulated  amission  from  an  uncorrelated  initial 
excited  state  in  a  gas  was  developed  through  a  calculation  of  the  transition 
probability  to  the  ground  state.  Further  complexity  was  introduced  into  the 
final  arrangement  of  photon,  thermal,  and  molecular  parameters  by  the  addition 
of  a  sound  wave.  In  particular,  the  precise  detuning  needed  to  produce  a 
maximuei  reeponse  of  stimulated  emission  directly  reflects  the  impact  of  the 
added  sound  wave  on  the  radiation  process.  Also,  a  thrse-level  model  was 
discussed  with  the  procedures  for  consideration  of  such  related  processes  as 
stimulated  Rasuin  scattering^. 

2.  To  explore  the  potential  of  a  traveling  current-initiated  (in  the  con¬ 
tainer  walls)  induced-dipole  laser  in  sMlecular  hydrogen  gas,  where  the  expec¬ 
tation  value  of  the  vector  potential  K  is  always  zero,  the  (purely  quantum; 
level  of  approximation  to  the  MECF  formalism  was  partially  developed.  It 
yielded  e  rigorous  fully  dimensional  traveling  wave  description  for  both  spon¬ 
taneous  emission  and  induced-dipole  emission  in  a  gas.  Specifically,  the 
ratio  of  the  rate  of  population  change  at  time  t  in  the  induced  dipole  case 

to  that  in  the  case  where  a  direct  non-induced  transition  can  occur  was 
directly  exhibited.  In  this  purely  quantum  example,  ^  can  no  longer  charac¬ 
terize  the  electromagnetic  field.  It  is  ectually  the  source  correlation  func¬ 
tions  which  couple  the  sMlecule  to  A  and  its  partial  derivative  that  determine 
the  Poynting  vector  and  so  characteFize  the  field.  (Results  have  not  been 
published.) 

3.  To  review  the  process  of  excitation  of  a  molecular  system  (hydrogen) 
by  electron  impect  (as  in  a  discharge  or  e-beam),  the  MECF  formalism  was 
displayed  for  these  two  mutually  interacting  species.  The  impact  of  pressure 
broadening  was  briefly  considered  and  a  group  of  new  matrix  elements  specifi¬ 
cally  involved  in  the  excitation  process  were  calculated  for  a  subset  of  the 
quantum  states  of  anlecular  hydrogen.  (Results  have  not  been  published.; 

Because  of  its  inherent  ability  to  separate  internal  structure  (symbolized 
by  the  wave  function  u)^  or  just  k)  from  center  of  mass  motion  (symbolized  by 
the  coordinate  x),  the  MECF  approach  is  particularly  applicable  to  processes 
such  as  RasMn  and  Brillouin  scattering  which  exhibit  contrasting  profiles 
relative  to  these  specific  properties.  Another  area  of  direct  potential  rele¬ 
vance  for  the  MECF  approach  is  the  investigation  of  Rydberg-state  interac¬ 
tions.  Here,  one  has  the  specific  ability  to  handle  a  complex  gaseous  mix 
involving  low  lying  states,  highly  excited  states,  ions,  electrons,  a  buffer 
gaa,  and  (of  course)  the  electromagnetic  field.  Also,  one  can  present 
numerous  other  applications,  both  standard  and  more  exotic,  such  as  the  impact 
of  light-  induced  drift  and  light  pressure  sound  waves  in  s  gas,  etc.  The  spe¬ 
cific  introduction  of  long-range  interaction  character,  introduced  through 
taking  the  transverse  part  of  S^j(x),  with  application  to  individual  atoms  and 
molecules  would  likewise  represent  a  potential  area  of  associated  interest. 


•Charl  as  Wang,  Ed.,  "Gas  Effects  in  Stimulated  Emission,"  in  Lasers  in 
Fluid  Mechanics  and  Plasmadynamics , "  AIAA,  New  York,  1983. 
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Finallyy  «  f«w  pot«nti«l  areas  of  application  that  axtand  the  limits  of 
the  formalism  are  mentioned.  One,  previously  alluded  to,  is  an  in-depth  con¬ 
sideration  of  the  interaction  of  electromagnetic  waves  with  a  fog  or  mist. 
Here,  one  has  increased  the  fundamental  sise  of  the  working  unit  from  that  of 
a  molecule  to  that  of  a  group  of  molecules,  etc.,  that  can  be  described 
collectively  as  having  both  a  center  of  mass  and  an  internal  character  or 
structure  analogous  to  that  of  a  large  swlecule.  Increasing  gaseous  pressure 
and  density,  or  the  other  hand,  eventually  smves  the  HECF  approach  toward  s 
consideration  of  transitions  in  phase.  Application  here  becomes  limited  when 
it  is  no  longer  possible  to  separate  out  the  individual  molecule.  Applica¬ 
tions  involving  a  consideration  of  lower-order  relativistic  effects  or  correc¬ 
tions  appear  to  be  sure  attractive  where  one  is  not,  as  yet,  forced  to  suve 
from  an  initial  Hamiltonian  form  to  a  consideration  of  alternative  Lagrsngisn 
formulations.  A  major  application  here  would  be  a  rigorous  characterization 
and  treatment  of  the  free-electron  laser,  where  the  basic  unit  is  now  the 
electron  and  relativistic  effects  have  been  introduced.  Surface  shielding 
effects,  involving  the  interaction  of  plasmas  with  radiation  or  particle 
beams,  represent  another  area  of  advanced  activity  treatable  by  the  MECF 
approach.  Here  the  specification  of  adequate  boundary  conditions  would  be 
particularly  challenging.  Additional  applications  can  also  be  cited  in  the 
general  area  of  collective  processes  and  phenomena. 
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APPENDIX 


An  issue  that  cannot  be  ignored  is  the  impact  of  the  Coulomb  interaction 
in  modifying  the  diagonal  elements  of  the  particle  density  matrix  in  the 
absense  of  any  external  or  applied  electromagnetic  field.  Any  modification  of 
this  nature  tiould  become  directly  applicable  to  our  study,  once  an  external 
field  was  added  to  the  system.  Because  this  modification  first  arises  at 
order  N^,  in  a  perturbation  expansion  on  the  Coulomb  interaction,  it  was  not 
directly  addressed  in  the  text,  which  chose  to  consider  only  terms  involving 
N^  and  N^.  However,  since  the  power  of  N  is  not  an  absolute  guide  to  (or 
measure  of)  impact  on  specific  physical  properties  and  observables,  an  analy¬ 
sis  that  can  be  used  toward  the  clarification  of  this  important  subject  is 
developed  in  this  Appendix. 

As  a  physical  example,  the  hygrogen  atom  is  again  chosen,  along  with  Its 
ground  state  1,  the  impact  of  the  three  levels  2(j)  that  produce  the  Lyman-a 
line  are  Investigated.  For  convenience,  at  this  time,  off-diagonal  correla¬ 
tion  functions  are  expressed  in  the  form  Rij(r^,t)  ■  Pl2(j)^£.»^^  1.  ^2(  j )  1  ^£.*  ^  ^  * 
which  due  to  the  absense  of  an  external  field  should  be  independent  of  the 
global  variable  Furthermore,  if  the  diagonal  elements  when  they 

appear  on  the  right  hand  side  of  the  equations  of  motion  are  allowed  to  assume 
the  seroth  order  form  gg,?^ O ,  or  a  simple  Gaussian  in  r,  and  then  transformed 
to  a  reduced  dimensionless  representation  (see  expressions  (16)),  the 
equations  of  motion  for  Fl.j(£,t)  are  simply 

i3tP'j(£,t)  -  Qi2P*j(lft)  ■  0, 


iatP^j(£,t)  -  ■  (l)2ii3/2g.2ni2Nj*d3kvj2(k)Q(-2-l/2^;2“l/2k)^ 

(A-1) 

where  the  dimensionless  (]i)  is  just  (k)  -  (g,ftn  V) 'l^g„,V^2(  j  H 
with  the  summation  over  ail  four  states, ^and  2a( E2( - )-Ej )7fi, 
g. ■gl-^2(j)  vhere  there  is  no  direct  dependence  on  j .  In  arriving  at  this 
simple  format,  the  fact  that  the  integral  over  all  space  of  Vg,2(  j  )ml 
sero  has  been  used  and  higher  order  terms  have  been  neglected  on  the  right 
hand  side  of  the  equations.  Since  time  t  and  ni2>  in  a  sense,  act  indepen¬ 
dently  in  this  development  where  they  are  directly  coupled,  they  may  be  con¬ 
sidered  to  be  in  either  the  original  or  the  reduced  dimensionless  form.  Also, 
the  function  Q,  as  it  appears  in  expressions  (A-1),  is  just  the  Q  of 
expression  ( 19) . 


As  a  measure  of  the  correction  on  the  diagonal  elements  Pmm(f>()  order 
N^,  the  quantity  AP(r_,t)  ■  APii(r^,t)  -  ?6P2(  j )  2(  j )  (£.»  t used,  which  folds 

all  of  the  three  2(j)  levels  into  the  mix.  In  this  way,  one  is  able  co 
illustrate  the  more  general  impact  of  the  Coulomb  interaction,  not  being  tied 
to  a  single  2(j)  level.  Using  the  same  approach  as  was  taken  in  the  construc¬ 
tion  of  expressions  (A-l),  one  can  now  write  the  equation  of  motion  for 
AP(r,t)  as 


i3cdP( 


t)  -  (  ♦)2g_^12Ny*d3kyd3x  exp(  ik'i^)  P(  x.)  ^  vj  ^  (ji)  j  (  x^>£,  t ) , 

(A-2) 


A-1 


where  only  P^j  is  directly  Involved.  Assuming  AP(£,0)"0,  expressions  (A-1) 
end  (A**2)  yield  a  simple  solution  with  oscillatory  time  dependence 

AP(£,t)  - 

P(r’)  2  7^2(31/2^)  7^2(31/2^') 

•  exp  t -(k2-k ‘k  •  +k •  2 )  ] exp  [ -1  (k-^' )  •  r  •  ] ,  ( A-3 ) 


where  £''"3''^/2^  end  P(£)  is  the  initial  Gaussian  functional  form.  Although 
expression  (A-T)  illusFrates  the  analytic  nature  of  AP(£it),  and  thus  and 

AP2(J)2(J)  through  the  obvious  relations  aPii*1/2aP,  JAP2(J)2(J)"  only 

by  inserting  the  latter  in  place  of  P(£)  in  expressions  such  as  the  trans* 
formed  version  of  (13)  can  one  truly  evaluate  its  quantitative  significance. 

It  is  obvious  from  (A-*!)  that  AP(£tt)  is  of  even  parity  in  r  and  is  saved  from 
approaching  sero  as  r  approaches  sero  by  the  cross  term  k*k^  that  couples  the 
k  and  k’  Integrals  together.  Further  details  on  this  subject  will  be  devel' 
oped  at  a  later  tine. 
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